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Abstract. We present observational estimates of the magnitude difference between the luminosity function red 
giant branch bump and the horizontal branch (AF555W5l^'^^''), and of star counts in the bump region (-Rbump), 
for a sample of 54 Galactic globular clusters observed by the HST. The large sample of stars resolved in each 
cluster, and the high photometric accuracy of the data allowed us to detect the bump also in a number of metal 
poor clusters. 

To reduce the photometric uncertainties, empirical values are compared with theoretical predictions obtained 
from a set of updated canonical stellar evolution models which have been transformed directly into the HST 
flight system. We found an overall qualitative agreement between theory and observations. Quantitative estimates 
of the confidence level are hampered by current uncertainties on the globular cluster metallicity scale, and by 
the strong dependence of AF555W^g™'' on the cluster metallicity. In case of the iibump parameter, which is 
only weakly affected by the metallicity, we find a very good quantitative agreement between theoretical canonical 
models and observations. For our full cluster sample the average difference between predicted and observed i?bump 
values is practically negligible, and ranges from —0.002 to —0.028, depending on the employed metallicity scale. 
The observed dispersion around these values is entirely consistent with the observational errors on -Rbump. As a 
comparison, the value of -Rbump predicted by theory in case of spurious bump detections due to Poisson noise in 
the stellar counts would be ~0.10 smaller than the observed ones. 

We have also tested the influence on the predicted AF555Wg™'' and -Rbump values of an He-enriched component 
in the cluster stellar population, as recently suggested by D'Antona et al. II2002II . We find that, under reasonable 
assumptions concerning the size of this He-enriched population and the degree of enrichment, the predicted 
AF555W^™^ and -Rbump values are only marginally affected. 
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1. Introduction 

Stellar evolution models supply fundamental tools which 
enable one to constrain the structure of the Milky Way, the 
properties of extragalactic stellar systems, and the early 
evolution of the universe. During the last few years the 
substantial increase in the spatial resolution provided by 
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the Hubble Space Telescope (HST), as well as the advent 
of wide field imagers in ground-based telescopes provided 
homogeneous and accurate photometry for large samples 
of stars in Galactic Globular Clusters (GGCs). This has 
made possible a thorough comparison between theoretical 
models of low mass stars and observations, over a broad 
metallicity range. This kind of comparisons play a crucial 
role in modern stellar astrophysics, because they provide 
stringent tests for the accuracy of current stellar evolution 
models (Renzini & Fusi Pecci .1988.: Castellani 1999). 

In this context, the Red Giant Branch (RGB) lumi- 
nosity function (LF) of GGCs is an important tool to 
test the chemical stratification inside the stellar envelopes 
(Renzini & Fusi Pecci il988,l . The most interesting feature 
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of the RGB LF is the occurrence of a local maximum in 
the luminosity distribution of RGB stars, which appears 
as a bump in the differential LF, and as a change in the 
slope of the cumulative LF. According to Thomas (1967) 
and Iben H1968fl this feature is caused by the sudden in- 
crease of H-abundance left over by the surface convection 
upon reaching its maximum inward extension at the base 
of the RGB. When the advancing H-burning shell encoun- 
ters this discontinuity, its efficiency is affected (sudden in- 
crease of the available fuel), causing a temporary drop of 
the surface luminosity. After some time the thermal equi- 
librium is restored and the surface luminosity starts to 
increase again. As a consequence, the stars cross the same 
luminosity interval three times, and this occurrence shows 
up as a characteristic peak in the differential LF of RGB 
stars. Moreover, since the H-profile before and after the 
discontinuity is different, the rate of advance of the H- 
burning shell changes when the discontinuity is crossed, 
thus causing a change in the slope of the cumulative LF. 

The brightness of the RGB bump is therefore related 
to the location of this H-abundance discontinuity, in the 
sense that the deeper the chemical discontinuity is located, 
the fainter is the bump luminosity. A comparison between 
the predicted bump luminosity and the observations al- 
lows a direct check of how well theoretical models for RGB 
stars predict the extension of convective regions in the 
stellar envelope; following Fusi Feed et al. ((1990 ) . the ob- 
served magnitude difference between the RGB bump and 
the HB at the RR Lyrae instability strip (AVj^™'^) is usu- 
ally employed in order to test the theoretical predictions 
for the bump brightness. This quantity presents several 
advantages from the observational point of view (see Fusi 
Pecci et al. [1990! and Salaris et al. I2U(J2|I and it is empiri- 
cally well-defined because it does not depend on a previous 
knowledge of the cluster distance. 

Dating back to its first detection in the LF of the GGC 
47 Tuc by King, Da Costa & Demarque (1985,) . the RGB 
bump has been the subject of several theoretical and ob- 
servational investigations (Fusi Pecci et al. 119901 Cassisi & 
Salarisl997, Alves & Saraiedini lTMl Ferraro et al. .1999. 
Zoccah et aLUnni Bergbusch & Vandenberg lMTTl Salaris, 
Cassisi & Weiss EOISI and references therein). However, a 
sound quantitative comparison between theory and obser- 
vations was hampered by the size of the observed stellar 
samples along the RGB, and by the heterogeneity of the 
datasets. This problem was even more severe for the most 
metal-poor GGCs, since the RGB evolutionary timescales 
is significantly shorter in metal-poor than in metal-rich 
stars. However, Zoccali et al. (1999, hereinafter Z99), us- 
ing a homogeneous set of data collected with HST, firmly 
detected the RGB bump in a large sample of GGCs cov- 
ering a wide metallicity range. 

This is the sixth in a series of papers aimed at in- 
vestigating the RGB bump in the LF of GGCs. On the 
basis of a detailed set of evolutionary models Cassisi & 
Salaris H1997|l showed that the predicted AV^g^'^was in 
agreement with empirical estimates of GGCs with accu- 
rate spectroscopic measurements of heavy-element abun- 



dances. This finding was further strengthened by the ev- 
idence that the AFj^g^^ is only marginally affected by 
atomic diffusion (Cassisi, DeglTnnocenti & Salaris [T997|l . 
A more detailed comparison was provided by Z99 who 
found that theory and observations do agree at the level 
of « 0.1 mag. 

A different kind of analysis was performed by Bono 
et al. 1)20011 hereinafter BOl), who compared the theoret- 
ical evolutionary lifetimes during the crossing of the H- 
discontinuity with the star counts across the RGB bump, 
a quantity that is sensitive to possible extra-mixing pro- 
cesses below the formal convective boundary. BOl have 
shown that star counts in the bump region are sensitive 
to the size of the jump in the H-profile left over when the 
envelope starts to recede, after achieving its maximum 
inward extension. Deep mixing phenomena (prior to the 
bump stage) able to dredge up an appreciable amount of 
He would change the size of the H-abundance jump, thus 
modifying the star counts in the bump region. The i?bump 
parameter was defined as the ratio between the star counts 
in the bump region Vbump ±0.4 and the star counts in the 
normalization region Vbump + 0.5 < Vbump < Vbump + 1-5. 
The reasons for this choice of the bump and the normal- 
ization regions are that the former should be large enough 
to include all bump stars, while the latter is required 
to normalize the total number of stars sampled in each 
cluster. The comparison between theory and observations 
performed by BOl disclosed that the occurrence of a sub- 
stantial efficiency of non-canonical extra-mixing before the 
RGB bump could be excluded. On the basis of theoretical 
arguments, Cassisi, Salaris & Bono 1)2002)1 more recently 
suggested to use the shape of the LF bump as a com- 
plementary diagnostic of partial mixing processes at the 
base of the outer convective zone, since it is sensitive to 
the shape of the H-discontinuity. 

The key differences between this investigation and Z99 
are the following: we adopted the sample of GGCs pre- 
sented by Piotto et al. 1)2002)1 and for 54 of them we de- 
tected the RGB bump. This sample is approximately 30% 
larger than the sample by Z99; ii) we adopted the more 
robust HST flight photometric system (without reddening 
corrections) instead of the Johnson de-reddened bands. 
This approach avoids deceptive errors in the estimate of 
the visual magnitudes of both RGB bump and HB. As a 
matter of fact, the calibration to the standard Johnson 
system fDolDhin .2000jl requires the knowledge of the red- 
dening. Therefore, the accuracy of the equivalent Johnson 
magnitudes depends on the accuracy of the adopted clus- 
ter reddening (see Piotto et al. l2002jl : in) the uncertainties 
of current parameters have been estimated using the com- 
pleteness curves and the photometric errors obtained from 
a large number of artificial star experiments (see iJ^J; iv) 
the comparison between theory and observations is based 
on an updated theoretical framework. 

The main aim of this investigation is to provide new 
homogeneous measurements of the RGB bump, and in 
turn of its brightness with respect to the HB, for a sizable 
sample of GGCs that cover a wide range in metallicity 
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(—2.2 < [Fe/H] < —0.3), and to compare these measure- 
ments with predictions based on a new set of evolutionary 
models. 

In the next section we present the observational 
database, and we briefly discuss the procedures adopted 
for the reduction and the calibration of the data. In this 
section we also outline the approach adopted to estimate 
the AF555WJ?[b"^ and the i?bump parameters. In §3, we 
present the new theoretical models adopted in this inves- 
tigation, while in §4 we compare theory and observations. 
Conclusions and future developments are discussed in §5. 

2. The cluster database 

We exploited our large photometric database of 74 GGCs 
observed in the HST B (F439W) and V (F555W) bands 
with the WFPC2 (Piotto et al. to measure the 

AF555Wy™^ and i?bump parameters, that are related to 
the position and extension of the RGB bump. The obser- 
vations, pre-processing, photometric reduction, and cali- 
bration of the instrumental magnitudes to the HST flight 
system as well as the artificial star experiments performed 
to derive star count completeness are described in Piotto 
et al. ((2D02). All photometric data have been processed 
following the same steps: after the pre-processing, the in- 
strumental photometry for each cluster was obtained with 
DAOPHOT II/ALLFRAME; the correction for the CTE 
effect and the calibration to the flight system was per- 
formed following the prescriptions by Dolphin H2()()()|l . 

For each cluster we measured a number of stars that 
ranges from a few thousands in the less massive clusters, 
to « 47,000 in NGC 6388. This and the high internal 
photometric accuracy, (the typical photometric error at 
the bump and ZAHB level ranges from a few 10"'^ to a 
few 10~^ mag) allowed us to detect the RGB bump even 
in metal-poor clusters. 

2.1. Measurement of the AFBbBW^^"^ parameter 

The parameter AFSSSWj?;'^^ used in this work is defined 
as the difference between the RGB bump and Zero Age 
HB (ZAHB) HST F555W magnitudes; it is therefore nec- 
essary to detect the bump along the cluster RGB and eval- 
uate the ZAHB brightness at the level of the RR Lyrae 
instability strip, following Z99. 

To measure the position of the RGB bump we used 
both the cumulative and differential LFs corrected us- 
ing the completeness functions for the RGB sequences 
obtained from the artificial star experiments (Piotto et 
al. 12002(1 . For those clusters from which more populous 
samples of RGB stars were obtained (e.g. the 26 clusters 
marked with an asterisk in Tab. 1) the cumulative LF 
alone provides a solid determination of the bump level, 
which is always consistent with the existence of a statis- 
tically significant bump in the differential LF (see Fig. 1 
for three examples). In case of the less populated clusters, 
both cumulative and differential LF have to be used si- 



multaneously in order to estimate the bump level, as first 
discussed by Fusi Pecci et al. (1990). As discussed in sec- 
tions 4.1 and 4.2, the bump detections in the more poorly 
sampled clusters are genuine, and in general, our detec- 
tions for the whole sample are internally consistent. 

We detected the RGB bump in 54 clusters, whereas for 
20 other clusters this was not possible due either to poor 
RGB statistics, or to a large differential reddening, or to 
high field star contamination, especially for bulge clusters. 

Fig. n shows some examples of the bump identi- 
fication in the differential and in the cumulative LFs 
for three clusters in our database: a metal-poor clus- 
ter (NGC 7078, [M/H]cG97 = -1-91), an intermediate- 
metallicity (NGC 6229, [M/H]cg97 = -1.09), and a 
metal-rich one (47 Tuc, [M/H]cg97 = -0.56). In aU pan- 
els, the position of the bump is marked by a vertical arrow. 

The error bars on the bump position have been esti- 
mated by adding in quadrature the photometric error at 
the bump magnitude and the half width of the bin size 
used in the RGB LF. 

Estimation of the ZAHB luminosity is a risky proce- 
dure due to the variety of HB morphologies and also to 
the fact that we cannot directly evaluate the mean magni- 
tude of the RR Lyrae stars from our photometry. In fact, 
our photometric data cover only a very short time interval 
and this does not allow a suitable sampling of the pulsa- 
tion cycle of the variables. Therefore, we followed the same 
procedure adopted by Z99. It can be brieffy summarized 
as follows: 

— clusters with [Fe/H] < —1 (metal- intermediate and 
metal-poor): we divided the sample into metallicity 
groups (each group spans at most 0.2 dex in metal- 
licity); to each metallicity group we associated a tem- 
plate cluster (with a metallicity within the group) ; the 
template clusters have accurate photometry and, most 
importantly, a sizable number of RR Lyrae, so that 
we could obtain an accurate estimate of the RR Lyrae 
magnitude level. Then we registered the HB of each 
cluster in the group to the HB of the reference clus- 
ter by artificially shifting the CMD of the template 
cluster in both color and magnitude in order to over- 
lap both their RGBs and HBs. The mean magnitude 
of the RR Lyrae for each cluster was then obtained 
by shifting the RR Lyrae magnitude of the reference 
cluster by the magnitude difference adopted to overlap 
the two CMDs. Finally, we transformed the RR Lyrae 
mean magnitude into the ZAHB magnitude following 
Cassisi & Salaris (1997; see Z99 and Recio-Blanco et 
al. "2003" for further details); 

— clusters with [Fc/H] > — 1 (metal- rich) : we first es- 
timated the photometric error (crpsssw) at the level 
of the HB using our set of artificial star experiments; 
then we determined the lower envelope of the HB stel- 
lar distribution. The ZAHB magnitude was eventually 
fixed at 3crF555W magnitudes above the lower envelope. 

The individual ZAHB luminosities as well as a more 
detailed description of the adopted procedure will be dis- 
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Table 1. Data for the 54 GGCs in our sample. The 26 clusters with more than 85 stars in the RGB Bump area are 
marked with an asterisk. 
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Fig. 1. Example of the bump detection in three clusters. Panel A - The two panels show the cumulative (top) and 
differential (bottom) RGB LF for the metal poor cluster NGC 7078; Panel 5 - As in Panel A, but for the intermediate 
metallicity cluster NGC 6229; Panel C - As in Panel A, but for the metal rich cluster 47 Tuc. The global metallicity 
on the CG97 scale (see the text for details) is shown in each panel. The bump position (marked by the arrow) clearly 
shows up as a change in the slope of the cumulative LF, and it is confirmed by the differential LF. 



cussed in a companion paper (Recio-Blanco et al. 2003, in 
preparation) . 

The errors in AF555Wj?[™^ have been estimated by 
summing in quadrature the errors in the bump position 
F555Wbump and in the ZAHB level F555Wzahb- 

Because of the uncertainty that still affects the metal- 
licity scale for GGCs (Rutledge, Hesser & Stetson 1997i 
VandenBerg HUM Caputo & Cassisi HUU^ and Kraft & 
Ivans 12003) I we have used both the Carretta & Gratton 
l(TM71 hereinafter CG97) and the Zinn & West 1(1^51 
hereinafter ZW84) metallicity scale. Moreover, we adopt, 
due to the lack of individual spectroscopic measurements 
for several GGCs in our sample, a mean a-enhancement 
of 0.3 dex for metal-poor and metal-intermediate clus- 
ters ([Fe/H] < —1.0) and of 0.2 dex for metal-rich clus- 
ters ([Fe/H] > —1.0). The former value was suggested 
by Carney ifT^ C96), while the latter is a mean be- 
tween the estimates collected by C96 and those collected 
by Salaris & Cassisi (1996). For each metallicity scale, the 
global cluster metallicity was estimated by using the re- 
lation provided by Salaris, Chieffi & Straniero Ijl993|l . We 
adopt a global metallicity error of ± 0.15 dex. This value 
can be considered as a safe lower limit to the uncertainties 
affecting both [Fe/H] and [a/Fe] measurements (Rutledge 
et alUMZIl- 

Table ^ lists the relevant data for all clusters in our 
sample where the RGB bump has been detected. Column 
(1) gives the identification, columns (2) and (3) give the 
cluster global metallicity according to the metallicity scale 
provided by CG97, and by ZW84, respectively. Column (4) 
lists the F555W magnitude of the RGB LF bump, column 
(5) lists the AFSSSW^"^^"^ values with associated errors. 



2.2. Measurement of the iibump paramenter 

We measured the i?bump parameter (using the same def- 
inition of BOl but the F555W magnitude) for all the 54 
GGCs where we detected the bump. The star counts in 
the bump and in the normalization regions were corrected 
using the completeness functions given by the artificial 
star experiments. 

The stellar counts in Bump and Normalization area 
are listed in columns (6) and (7) of Tab. 1, respectively. 
The values of the i?bump parameter with the correspond- 
ing errors are listed in column (8). The la error on 
the observed i?bump values is determined according to 

CT(i?bump)=-Rbump Vi^/^b) + (l/iV„) whcrc Nb and iV„ 
are the star counts in the bump and in the normalization 
region, respectively. 

3. The theoretical framework 

The theoretical framework we adopt in this investigation 
is based on an updated and larger set of stellar models 
(Pietrinferni, Cassisi, Salaris & Castelli 2003). The new 
models have been computed by using a recent version of 
the FRANEC evolutionary code. The input physics has 
been updated with respect to the models used in Z99 
(Cassisi & Salaris 1997), and the changes are summarized 
in the following: 

— the radiative opacity is obtained from the OPAL tables 
(Iglesias & Rogers il996i for temperatures larger than 
10'' K, and from Alexander & Ferguson l|1994|l for lower 
temperatures. Opacity for electron degenerate matter 
is computed following Potekhin l|1999|l . 

— We updated the energy loss rates for plasma-neutrino 
processes by using the most recent and accurate re- 
sults provided by Haft, Raffelt & Weiss fjl994)l . For all 
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other processes we still rely on the same prescriptions 
adopted by Cassisi & Salaris (1997). 

— The nuclear reaction rates have been updated by using 
the NACRE database (Angulo et al. 1999), with the 
exception of the ^^C(a,7)^^0 reaction. For this reac- 
tion we employ the more accurate recent determination 
by Kunz et al. PUJ^ . 

— The detailed Equation of State (EOS) by A. Irwin^ 
has been used. A full description of this EOS is still 
in preparation (Irwin et al. 2003) but a brief discus- 
sion of its main characteristics can be found in Cassisi, 
Salaris & Irwin (2003). It is enough to mention here 
that this EOS, whose accuracy and reliability is sim- 
ilar to the OPAL EOS developed at the Livermore 
Laboratories (Rogers, Swenson & Iglesias '19961 and 
recently updated in the treatment of some physical in- 
puts (Rogers & Navfonov l2002|l . allows us to compute 
self-consistent stellar models in all evolutionary phases 
relevant to this investigation. 

— The extension of the convective zones is fixed by means 
of the classical Schwarzschild criterion. Induced over- 
shooting and semiconvection during the He-central 
burning phase are accounted for following Castellani et 
al. H1985|l . The thermal gradient in the superadiabatic 
regions is determined according to the mixing length 
theory, whose free parameter has been calibrated by 
computing a solar standard model. 

— Our evolutionary code includes the process of atomic 
diffusion of both Helium and heavy elements; the mod- 
els used in this work, however, have been computed 
with the atomic diffusion switched off, taking into ac- 
count that Cassisi et al. H1997|l have clearly shown how 
the effect of this physical process on Ayjjg™''is negli- 
gible. 

— We adopt a scaled-solar heavy element mixture 
(Grevesse & Noels '1993|; the enhancement of a- 
elements observed in GGCs is accounted for following 
the prescriptions by Salaris et al. I|1993p . 

— As far as the initial He-abundance is concerned, we 
adopt the estimate recently provided by Cassisi et al. 
120031 on the basis of new measurements of the R pa- 
rameter in a large sample of GGCs. They found an 
initial He-abundance for GGC stars of the order of 
Y = 0.245, which is in fair agreement with recent em- 
pirical measurements of the cosmological baryon den- 
sity provided by W-MAP (Spergel et al. 2003). To re- 
produce the calibrated initial solar He-abundance we 
used dY/dZ w 1.4 (Pietrinferni et al. 2003). 

— Bolometric magnitudes have been transformed into 
HST F555W magnitudes according to the transforma- 
tions provided by Origlia & Leitherer l)2000|l , based on 
the atmosphere models computed by Bessell, Castelli 

& Plez (cnnBj- 



^ The EOS code is made publicly available at 
fftp: / / astroftp.phy s.uvic.ca| under the GNU General Public 
License (GPL) 



Table 2. Theoretical AF555Wh™p values as a function 
of [M/H] and age. 



[M/H] AF555W^™P 





10 Gyr 


12 Gyr 


14 Gyr 


16 Gyr. 


-2.267 


-0.907 


-0.846 


-0.770 


-0.699 


-1.790 


-0.690 


-0.621 


-0.555 


-0.494 


-1.266 


-0.297 


-0.221 


-0.159 


-0.116 


-0.963 


-0.056 


0.004 


0.056 


0.108 


-0.659 


0.192 


0.271 


0.347 


0.418 


-0.253 


0.469 


0.566 


0.649 


0.694 


0.000 


0.639 


0.741 


0.791 


0.843 



A more detailed discussion about the model input physics 
and the evolutionary code is presented in Pietrinferni et 
al. (2003). 

Finally, we discuss in more detail the metal distribu- 
tion used to compute the evolutionary models. In our 
analysis we will compare the models computed with a 
scaled-solar metal abundance to GGCs data for which 
the metal abundance is a-enhanced. According to Salaris 
et al. (1993), as long as the a-elements have approxi- 
mately the same enhancement, scaled-solar models mimic 
a-enhanced models computed with the same global metal- 
licity [M/H]. This is, however, strictly true only for 
[M/H] values up to ^ —1.0, as shown by Salaris & 
Weiss (1998) and Vandenberg et al. (2000). For higher 
global metallicities this equivalence is not well satisfied 
anymore. Therefore we investigated how different would 
be the AFSSSWj^jg"'^ parameter at a given [M/H] larger 
than —1, for a scaled-solar and an a-enhanced metal 
mixture. We used isochrones by Salaris & Weiss (1998) 
with [M/H]=— 0.3 (approximately the upper end of the 
metallicity range spanned by our GGC sample), both 
scaled-solar and with [a/Fe]=0.4; we found that the 
AF555Wy™'' values are changed by only 0.05 mag for 
typical GGC ages, the a-enhanced ones being larger. This 
is however just an upper limit to the real difference, since 
these metal rich clusters seem to show an a-enhancement 
lower than [a/Fe]=0.4. We conclude that, as far as the 
AF555Wh™^ parameter is concerned, we can safely use 
scaled-solar models with the cluster global [M/H] even 
in the high metallicity regime. We also verified that the 
-Rbump values are unaffected by the a-enhanced metal dis- 
tribution. 

4. Comparison between theory and observations 

4.1. The AF555WhTp parameter 

Figure El shows the comparison between theoretical pre- 
dictions (given in Tab. 2) and measured values of 
AF555Wjj™'' as a function of the global metallicity 
[M/H], for a selected age range. The behaviour of the the- 
oretical AF555Wj?['^^'' stems from the strong dependence 
of the depth of the convective envelope on both metal- 
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Fig. 2. Comparison between the theoretical and empirical values of AF555Wy™^ as a function of the global metallicity 
for all the 54 GGCs in our sample (left panel) and for the 26 GGCs with more than 85 stars in the bump region (right 
panel). The empirical AF555Whb™p data are plotted as a function of both the Zinn & West l|1984l - ZW84) and 
Carretta & Gratton H1997I - CG97) metallicity scales (see text for more details). The theoretical predictions are 
plotted for four different cluster ages, namely 10 Gyr (solid line), 12 Gyr (long dash), 14 Gyr (short dash), and 16 Gyr 
(dot - short dash). 



licity and mass. The higher the cluster age, the smaller 
the evolving mass along the RGB, with an ensuing deeper 
extension of the convective envelope, and a fainter bump 
brightness, whereas the HB level is largely unaffected. On 
the other hand, at a fixed age, the larger the metallicity, 
the deeper the extension of the convective envelope which 
causes a fainter bump; the HB is also fainter for increasing 
metallicity, but the effect on the bump level is larger, thus 
causing an increase of AF555Wy™''. 

The left panel of Fig. 2 shows the data for the entire 
sample of GCs, whereas in the right panel we plot only the 
clusters with more than 85 stars in the bump region (here- 
after we will refer to this sample of clusters as the best 
cluster sample). A glance at this figure reveals that the 
distribution of the observational data for the entire sample 
closely resembles the distribution for the best cluster sam- 
ple which, incidentally, spans the entire metallicity range 
covered by the full sample. This evidence supports the sig- 
nificance of the bump detections for the whole 54 clusters 
listed in Tab. 1. Further confirmation comes from the fact 
that, as we will see in the following, the comparison with 
theory provides exactly the same results when considering 
either the full sample or the best cluster sample. We antic- 
ipate that the analysis of the i?bump parameter, presented 
in the next subsection, will provide additional evidence for 
the significance of the bump detection in all the clusters 
displayed in Table 1. 



Figure 13 shows also a qualitative agreement between 
observations and theory. In particular, both theoreti- 
cal lines and observations lie in the same region of the 
AF555W5i^™''[M/H] plane, and the observed trend with 
metallicity is well reproduced by theory (with the possi- 
ble exception of the three most metal poor clusters). 

We have also performed a more quantitative compari- 
son, though it must be noted that it is made difficult (and 
uncertain) by the strong dependence of AF555Wjjg™^ on 
the cluster age and metal content. Indeed, AF555Wyg"^^ 
varies by 0.03 magnitudes for a 1 Gyr variation in age, 
and by the same amount for a variation of only 0.04 dex 
in [M/H]. To assess the internal accuracy of theoretical 
models we compared the cluster age required to fit the ob- 
served AF555Wyg"'' values with the age determinations 
obtained from the main sequence (MS) turn off (TO) po- 
sition. It must be explicitly noted that we do not propose 
to use the AF555Wyg™'' parameter as an age indicator; 
however, we should expect that, for any given metallicity 
scale, the age scale inferred from the AF555W^™'' agrees 
with the ages from the TO. Systematic or random age 
offsets imply (systematic or random) effects not properly 
accounted for by the models. 

We started with the CG97 metallicity scale. The the- 
oretical line which best fits the observed distribution of 
AF555Wg™'' with the metallicity has an age of 11.8 Gyr, 
with a Icr dispersion of 4.0 Gyr. Figure^ (upper panel) dis- 
plays the comparison between the observed AF555W^™^ 



8 



M. Riello et al.: The Red Giant Branch Bump 



and the models for the labeled age. The use of only the 
best cluster sample (26 objects) sample does not affect 
the age, though the spread is slightly reduced. A close in- 
spection of Fig.|31shows that clusters with [M/H]< —1.6 
are systematically discrepant, even those with the best 
populated bump. This fact may point out to some serious 
inaccuracy of the stellar models for this metallicity range, 
though we cannot exclude that the discrepancy may be 
ascribed to a different 'true' global metallicity for these 
clusters (e.g., the a-enhancement is higher than assumed). 
Moreover, the discrepant clusters in this metallicity range 
are only a few, namely M 15, M 68, and M 53. A larger 
sample of metal-poor clusters with a well-populated RGB 
is required before a definitive conclusion on this possible 
discrepancy can be reached. 

If we restrict our analysis to the clusters with [M/H]> 
— 1.6, the best fitting model corresponds to an age of « 11 
Gyr, which does not change when we consider only objects 
belonging to the best cluster sample. The dispersion of the 
data points around this best fitting model does, however, 
decreases from 3.4 to 2.6 Gyr when using the best sam- 
ple. Most interestingly, we also found that the age scatter 
around the best fitting model does not correlate signifi- 
cantly with [M/H]. 

A different result is obtained using the ZW84 metal- 
licity scale. In this case, the age of the best fitting model 
is « 15 Gyr (Fig. O, but now the residuals show a cor- 
relation with [M/H]. In order to obtain a better match 
with the theoretical AF555Whb"''' 

a younger age for more 
metal-rich clusters must be used. More specifically, we 
find an overall statistically significant age- [M/H] relation- 
ship, with a slope (—5.7 ± l.l)Gyr dex^^. If we neglect 
the three most metal-poor clusters mentioned above (they 
have [M/H]< —1.8 in the ZW84 scale), we obtain a slope 
(— 4.8±l.l)Gyr dex^^, which does not change for the best 
cluster sample (with [M/H]> —1.8). In this case, the best 
fitting model correspond has an age of « 14 Gyr. 

In summary, when using the CG97 scale the cluster av- 
erage age is of 11-12 Gyr, in agreement with the TO ages 
obtained by Salaris & Weiss (2002), who determine the 
GGC distances from their theoretical models. Using the 
same metallicity scale, and measuring distances with the 
classical MS fitting technique, Carretta et al. (2000) find 
a similar age (though their best age estimate is roughly 13 
Gyr, when they consider also other independent distance 
scales) ^ . A second result of the comparison shown in Fig. 
3 is that, when using the CG97 metallicity scale, we do 
not find any statistically significant age-metallicity trend 
from the AF555Wji^™^ parameter, whereas Rosenberg et 
al. (il999), VandenBerg (2000), Salaris & Weiss (2002), 
find in general that the more metal rich clusters are 

^ Although a detailed discussion of the GGCs ages is clearly 
beyond the scope of this investigation, it has to be mentioned 
that a very recent estimate (Gratton et al. I2()(),'i|l of the abso- 
lute age of one metal-intermediate and one metal-poor cluster, 
based on local subdwarfs, and on new, accurate determinations 
of reddening and metallicity, is between 13.5 and 14 Gyr. 



younger, and Gratton et al. H20U3fl find that the metal- 
rich cluster 47 Tuc is about 2.6 Gyr younger than more 
metal-poor clusters like NGC 6752. 

Salaris & Weiss (2002) obtain, on average, ages higher 
by about 1 Gyr when moving from the CG97 to the ZW84 
metallicity scale, but still lower than the age derived from 
the AF555Wh™p. VandenBerg (2000) determines ages of 
about 14 Gyr for the oldest metal-poor GGCs using his 
own HB models and the ZW84 scale, and ages decreasing 
on average when moving to the higher metallicity regime. 
These results confirms the earlier work by Rosenberg et 
al. (19S9|, who find that the most metal rich clusters 
are 15-20% younger than the most metal poor ones. In 
conclusion, the age-metallicity relationship we found from 
AF555Wy™^ when adopting the ZW84 metallicity scale, 
is in qualitative agreement with the results we obtain from 
the TO ages, but the slope is larger. Also the average age 
implied by the AF555Whb"''' 

parameter seems larger than 
the most recent determinations. 

Unfortunately, this analysis provides somewhat not 
definitive results about the accuracy of the theoretical 
AF555Wy™'' values. This is due mainly to the current 
uncertainties in the cluster metallicities and, to a minor 
extent, to uncertainties in the GGC ages (which depend 
themselves on the metallicity errors). In fact, typical dif- 
ferences among the TO absolute GGC ages mentioned be- 
fore are of the order of 1-2 Gyr; these differences affect 
AF555Whb"''' (at a given [M/H]) at the level of 0.03-0.06 
mag, whereas typical differences of 0.2 dex between the 
CG97 and ZW84 [M/H] scale modify AF555WhT'' (at a 
fixed age) by ~0.2 mag. 

4.2. The i?bump parameter 

The comparison between theoretical predictions and em- 
pirical measurements of the i?bump parameter is shown in 
Fig-E for the same two metallicity scales discussed before. 
It is important to remark that, as in case of AF555Whb"'p 
the i?bump distribution for the clusters with the most pop- 
ulated bump region does not show any evident systematic 
offset when compared to the full sample. Due to the fact 
that one expects a weak dependence of i?bump on both 
[M/H] and age, it is meaningful to compare the mean 
measured i?bump values for the best cluster sample and 
the full sample. The mean value of i?bump for the latter is 
0.498±0.092 (Ict dispersion); in case of selecting only the 
26 objects belonging to the best cluster sample (which cor- 
responds to clusters with cr(i?bump)/-Rbumpbelow 13%) we 
obtain 0.497±0.064, practically identical to the value for 
the whole sample, but with a smaller dispersion because of 
smaller individual error bars (due to Poisson fluctuations 
in the values of Ni, and Nn). This is another clear evi- 
dence for the significance of our bump detections. In fact, 
in case of spurious detections in the less populated clus- 
ters, caused by random fluctuations in the number of RGB 
stars, we would expect to obtain an average i?bump value 
equal to ~0.37 (corresponding to models where the H- 
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[M/H] 

Fig. 3. Comparison between observed and predicted 
AF555Wh™p values for the labeled average ages deter- 
mined in case of, respectively, the CG97 (upper panel) 
and ZW84 (lower panel) [M/H] scale. Filled circles denote 
the 26 GGCs with more than 85 stars in the bump region, 
whereas open circles represent the remaining clusters in 
our sample. 




-2.0 -1.5 -1.0 -0.5 



[M/H] 

Fig. 4. Comparison between the theoretical and empirical 
values of i?bump as a function of the global metallicity for 
all the GGCs in our sample; symbols are as in Fig. |3| 
The bottom and upper panels refer, respectively, to the 
ZW84 and CG97 [M/H] scale. The sohd lines show the 
theoretical values for 10 (upper line) and 16 Gyr. 



abundance jump is absent) with a given dispersion around 
this value due to Poisson noise; however, when we consider 
only the less populated clusters with cr(i?bump)/^bump 
larger than 13%, we obtain a distribution of points with a 
mean value of 0.498 equal to the best cluster sample, and 
not centred around 0.37. 

To test the consistency of the observed i?bump with 
theory, we adopted the following approach: for a given 
metallicity scale we computed, on a cluster by cluster ba- 
sis, the corresponding theoretical i?bump values, using the 
ages obtained from AF555Whb"p. The precise age value 
is however not crucial, due to the very low sensitivity of 
the i?bump parameter to age (see the discussion in BOl and 
data plotted in Fig.^. Also the dependence on [M/H] is 
very weak, as noticed before. 

These 54 theoretical values have been compared 
with the empirical data and the difference between the- 
ory and observations has been analyzed (see Fig. [SJ. 
We found an average difference (theory-observations) 
A(i?bump)=-0.002±0.093 (Icr) for the whole GGC sample 
when using the CG97 metallicity scale, with no evident 
correlation of the residuals with [M/H]. We also tested, 
making use of Monte Carlo simulations, if the dispersion 
around the mean can be explained by the Poisson errors 
on the empirical determination of i?bump- More in detail, 
for each cluster we randomly generated 10000 values of 



A(i?bump), centred around —0.002 and with a Gaussian 
la dispersion equal to the observational error on i?bump- 
We then joined together the synthetic A(i?bump) values for 
each individual cluster, and determined the Icr dispersion 
of the resulting distribution, which compares well with the 
observed value of 0.093. We also considered, as a further 
test, the subsample of 26 clusters with the best populated 
RGB up to the bump region. The i?bump residuals show 
again a negligible difference from zero, the average value 
being — 0.002±0.065, and once again we found that the 
dispersion around this value can be explained as obser- 
vational errors on the i?bump parameter. Note the much 
smaller value of the dispersion when only the best i?bump 
determinations are accounted for. We obtain the same re- 
sult if M 15 and M 53 - the two most metal poor clusters 
- are neglected (see Fig. EJ. 

When using the ZW84 scale, and the age distribution 
obtained from AF555Wh™p, the average value of the dif- 
ference between theoretical and observed iibump is slightly 
higher, namely — 0.028±0.093, the residuals showing no 
trend with [M/H]. The average value is -0.022±0.064 if 
we consider the best cluster sample. 

4.3. The effect of chemical self pollution 

To explain the observed CNO abundance anomalies and 
the extended blue HB tail in some GGCs, D'Antona et 
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Fig. 5. Difference between tlie tlieoretical i?bump values 
computed for the individual cluster ages determined from 
AF555Wy™'' and the observational counterpart, in case 
of the CG97 (upper panel) and ZW84 (lower panel) [M/H] 
scale. The horizontal dashed line corresponds to a zero 
difference and it is not a fit to the data. 



al. (2002) suggested the existence of an He-enriched stel- 
lar component within the clusters, due to chemical pol- 
lution by the ejecta of massive asymptotic giant branch 
stars. We tested whether this He-enhancement might af- 
fect predicted AF555Wy™P and i?bump values when com- 
pared to the case of a constant He-abundance. Following 
D'Antona et al. (2002) results, we have produced a syn- 
thetic CMD for an hypothetical 12 Gyr old GGC with 
[M/H] =-1.3. We considered that 64% of the coeval stel- 
lar population is composed of stars with the standard He- 
abundance adopted in our models, whereas 36% of the 
population is made of stars with Y randomly distributed 
between the standard value and an abundance 0.06 higher. 
We found that the location of the bump in the LF of this 
He-enhanced population does not show any significant dif- 
ference with respect to the standard case. The HB level at 
the RR Lyrae instability strip - as shown by D'Antona et 
al. (2002) - is determined by the 'He-normal' population, 
and therefore, the AF555Wy™^ value for a cluster with 
He-enhanced stars is expected to be very close to the case 
of a standard He-normal population. It is also important 
to notice that, if the He-enhanced stars have all the same 
abundance 0.06 higher than the bulk of the cluster popu- 
lation, the cluster LF would show a second bump, about 
0.13 mag brighter than the main one. 

As far as the i?bump parameter is concerned, we found 
only a a marginal difference. In conclusion, the occurrence 



of chemical pollution due to the same GGC stars does not 
affect the AF555Wy™^ and the i?bump values predicted 
by canonical stellar models. 

5. Summary and final remarks 

In this investigation we adopted a database of HST pho- 
tometric data for a large sample of GGCs, and we have 
been able to measure the RGB bump location in a sample 
of 54 clusters. For the same cluster sample we have also 
determined the star counts in the bump region following 
BOl. This represents, so far, the largest database of em- 
pirical estimates for both the RGB bump brightness and 
star counts in this evolutionary phase. 

The observed magnitude difference between the bump 
and the HB, i.e. AF555Whb"p and the star counts in the 
bump region, i.e. the i?bump parameter, have been com- 
pared with theoretical predictions by using a new set of 
stellar evolution models. To account for the current un- 
certainty in the GGC metallicity scale, we have adopted 
both the CG97 [Fe/H] scale and the ZW84 one. We also 
employed reasonable assumptions for the a-element en- 
hancement in GGCs stars, based on the presently available 
estimates. 

Owing to the sensitivity of the AF555Wj?[g"'' param- 
eter on cluster age, the ages required to fit the observed 
AF555Wh™p values should agree with recent indepen- 
dent estimates based on the luminosity of the TO, for the 
theoretical models being consistent with observations. Our 
results can be summarized as follows: 

— The mean age needed to obtain agreement between 
observed and predicted AF555Wj?[™'' values is ~ 12 
Gyrs (with a Icr dispersion of 4.0 Gyrs) when the CG97 
[Fe/H] scale is used. This value is slightly smaller than 
recent independent estimates obtained by using the 
same metallicity scale. It is also worth mentioning that 
a significant discrepancy exists for the most metal-poor 
clusters ([M/H]< —1.6) in our sample. Current data do 
not allow us to conclusively assess whether this rep- 
resents a real problem for the theoretical models, or 
the discrepancy is due is an observational bias due to 
the limited sample of stars located along the RGB. 
At variance with results from studies of TO ages, the 
AF555Wh™p of metal rich clusters like 47 Tuc IS re- 
produced by models with the same average age as more 
metal poor ones. 

— The mean cluster age is ~ 15 Gyr when the ZW84 
[Fe/H] scale is used and, in addition, we find a sta- 
tistically significant correlation between the residuals 
around this age and the cluster [M/H]. This mean 
age of 15 Gyr appears larger than similar mean GGC 
ages available in the literature. The slope of the age- 
metallicity relationship, though in qualitative agree- 
ment with current estimates (e.g., Rosenberg et al. 
119991 and Gratton et al. 2003), is also too steep. 

On the basis of these results we can draw the follow- 
ing conclusion: even though a qualitative agreement be- 
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tween theory and observations of AFSSSWjj"^^^ does ex- 
ist, a more definitive assessment of the confidence level is 
hampered by the not negligible uncertainties still affecting 
both the cluster [Fe/H] and [a/Fe] estimates. More robust 
conclusions require more accurate spectroscopic measure- 
ments of these quantities. 

As far as the -Rbump parameter is concerned, our anal- 
ysis suggests that there is a good agreement between the- 
ory and observations, regardless of the adopted metallicity 
scale. This result is also due to the weak dependence of 
this parameter on the cluster metallicity (and age), which 
minimizes the effects related to the uncertainties on the 
[M/H] scale. 

Finally, we found that the effect on both AF555Wh™p 
and i?buinp of a possible He-rich component - as recently 
suggested by D'Antona et al. (2002) - in the cluster stellar 
population is negligible. 
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